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ABSTRACT 

C**- ' Astrometric and radial-velocity planet detections track very similar motions, 

and one generally expects that the statistical properties of the detections would 
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also be similar after they are scaled to the signal-to-noise ratio of the underlying 



observations. I show that this expectation is realized for periods small compared 
to the duration of the experiment P/T <C 1, but not when P/T > 1. At longer 
periods, the fact that models of astrometric observations must take account of 
an extra nuisance parameter causes the mass error to begin deteriorating at 
P/T ~ 0.8, as compared to P/T ~ 1.0 for RV. Moreover, the deterioration is 
much less graceful. This qualitative difference carries over to the more compli- 
cated case in which the planet is being monitored in the presence of a distant 



companion that generates an approximately uniform acceleration. The period 
errors begin deteriorating somewhat earlier in all cases, but the situation is quai- 
itatively similar to that of the mass errors. These results imply that to preserve 



astrometric discovery space at the longest accessible orbits (which nominally 
have the lowest-mass sensitivity) requires supplementary observations to identify 
or rule out distant companions that could contribute quasi-uniform acceleration. 

Subject headings: planetary systems - methods numerical 



1. Introduction 

Astrometric and radial-velocity (RV) planet detections are, from a mathematical stand- 
point, extremely similar. In each, one models 1-dimensional projections of Kepler orbits and 
attempts to fit Kepler parameters. In the limit of circular orbits, the planet signature is 
encoded in the simple form 



F(t; Oi, a 2 , a 3 ) = a 1 sin(a 2 t + a 3 ) 
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where 

a± — a ( astro metry), a\ — K (RV), (2) 

are the astrometric and velocity semi-amplitudes for the two cases, and where 

(i2 = — , 03 = (astrometry + RV) (3) 

designate the period P and phase 0, respectively. One then goes on to combine information 
about the star's mass and (for astrometry) its distance, to infer the planet mass m (astrom- 
etry) or msini (RV), where i is the inclination. In the astrometric case, there are of course 
two such equations, one for each direction in the plane of the sky, whose ratio gives cosz, and 
so permit one to break the msini degeneracy that plagues the intrinsically 1-dimensional 
RV measurement. 

Because the form of equation (pQ) is essentially identical in the two cases, it is generally 
assumed that the error properties, i.e., the relation between the measurement errors and the 
derived-parameter errors, is also the same. Of course, it is well known that the parameter 
errors have a different dependence on semimajor axis a, system distance D, etc. Most 
notably, with other parameters held fixed, astrometric sensitivity increases linearly with a 
whereas RV sensitivity declines as cT x l 2 . But here I am referring to something else. The 
differences just mentioned all impact the final result because they change the characteristic 
signal-to-noise ratio of the experiment, 

SNR=— y/N, (4) 

a 

where N is the number of measurements and a is the error in each measurement (assumed for 
simplicity to be all the same). Here I will show that the astrometric and RV measurements 
have substantially different error properties even when SNR is identical. 



2. Minimum Variance Bound 

I will work within the framework of the minimum variance bound (MVB), also frequently 
called the Fisher-matrix approximation. As the same approximation will be applied to both 
techniques, this will allow me to highlight the difference between them. Of course, if equation 
([I]) really did fully represent both techniques, there could not be any difference in their error 
properties. However, the true functional forms of the source motions are actually 

n 

F(t; di) = ai sin(a 2 t + a 3 ) + ^J a,^ 4-4 - (5) 

i=A 
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For simplicity, I will continue to assume circular orbits, and I will assume for the astrometric 
case that the orbit is seen edge-on. In fact, the latter is not much of an approximation since 
in most cases the great majority of the information about the mass and period comes from 
the major axis of the apparent astrometric ellipse. Finally, I will assume that the data are 
taken uniformly in time, at intervals that are frequent compared to the orbital period, and 
that they all have the same error, a. 

For a single-planet, i.e., a system composed of just a star and a planet, in which the star 
does not suffer any other accelerations, n = 4 for the RV case, but n = 5 for the astrometric 
case. For RV, a^ is the systemic radial velocity of the system. For astrometry, a^ is the 
systemic proper motion of the system, while a 4 is the zero point of that motion. 

The introduction of the nuisances parameters a 4 or (a 4 , a 5 ) is what distinguishes the two 
cases statistically. The difference appears only for P ~ T or P > T, where T is the duration 
of the experiment. I have implicitly assumed that T, P ^> 1 yr, and this is the reason that 
I ignore parallax, which would be an additional nuisance parameter for astrometry, but not 
for RV. 

In the Fisher-matrix approximation, the inverse covariance matrix of the parameters, 
kj, is given by 

in^E^P - " £ fdtfmm (6) 
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where 



One easily finds that for P ^ T, a\ and a 2 are completely uncorrelated from the other 
parameters, so that 
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<x\ a-y \fb~i~i~ai V N a\ SNR 

That is, a(ai)/ai, which is the fractional error in astrometric or RV amplitude depends only 
on SNR, and in a very simple way. Similarly, 

<r{P) _ ^M __ 1 = v^/tt P 
P a 2 ~* y/b^a 2 SNR T' [ } 

However, as P approaches T, the parameters of primary interest (mass and period) 
become correlated with the nuisance parameters. To calculate this effect, I employ equations 
(JSJ) and ([7j), but using the full 7-parameter Kepler formalism, not just the three Kepler 
parameters displayed in equations (TTT) and (EJ. That is, even though the adopted orbits 
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are circular and edge-on, I allow for free fits to the eccentricity and inclination, and so for 
correlations between these parameters (as well as the two remaining Kepler parameters) and 
the parameters of interest (amplitude and period) . The effect of allowing for these covariances 
is to increase the errors in semi-amplitude and period by modest amounts relative to what 
would be obtained using equation ([5j). 

Panels (a) and (b) of Figure [1] show the ratios of the true errors for the amplitude and 
period, respectively, relative to the naive equations (|Sj) and (|H]), for both the RV (green) 
and astrometric (red) cases. In fact, as I will discuss in § HJ for P > T, the errors depends 
on phase as well as period. Figure [1] therefore shows the root-mean-square of the errors, 
averaged over all phases. Note that the RV amplitude errors follow the naive form until 
P/T ~ 1.1 and then deteriorate relatively gracefully. By contrast, the astrometric errors 
begin deviating at P/T ~ 0.75 and then deteriorate much more quickly. For the period 
errors, deterioration begins at P/T ~ 0.85 for RV and P/T ~ 0.65 for astrometry, but the 
overall pattern is qualitatively similar. 

The mass estimates for astrometry and RV depend on different combinations of ampli- 
tude and period, 

/ 2 1 \ 

m oc aia'2 (n=- Ast; n = — RV] (10) 

Hence the fractional error in the mass (or m sin i in the case of RV) is related to the errors 
in the fit parameters by 

^ + 2n^ + n 2 ^. (11) 

In the limit P ^ T, a(m)/m — > a(a±)/ai, but for P > T, the period error and the corre- 
lations become important. Figure [lb shows the results of calculations that apply equation 

(hid. 
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3. Additional Uniform Acceleration 

Of course, the star may have more than one companion (planetary or otherwise), and 
one may imagine arbitrarily complicated configurations. Here I restrict myself to the next 
level of complication, a second companion that is sufficiently far away that its effect on the 
star may be treated as uniform acceleration. Even if no such acceleration is identified, one 
might decide to fit for it on the grounds that there may be such a companion that has not 
been recognized. I now ask how the inclusion of such an additional nuisance parameter 
affects the precision of the physical parameters for the planet that was previously treated as 
isolated. 
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The number of nuisance parameters is incremented by one in each case. Referring to 
equation (jSJ), for RV, we now have n — 5, with a§ being the systemic acceleration, while for 
astrometry, n — 6, with a 6 being the systemic proper-motion acceleration. Figures dh and 
[Tb show the corresponding ratios relative to the naive equations (jSj) and Q). In this case, 
RV is represented by the red curve, while astrometry is represented by the blue curve. The 
fractional mass errors are shown by solid curves in Figure [lb. 

In this case, both deteriorate rapidly, but the mass-error deterioration begins at P/T ~ 
0.65 for astrometry and at P/T ~ 0.75 for RV. One may anticipate that in other, yet more 
complicated situations, RV will always perform like the naive equations to higher P/T than 
astrometry, simply because it has one fewer nuisance parameter. 



4. Phase Dependence 

Figure [T] shows the root-mean-square errors averaged over phase. In fact, as P/T grows 
and the rms errors deteriorate, the variations with phase also increase. Figure [2] shows the 
phase variations of the mass errors for the particular case that the period is exactly equal to 
the duration of the experiment, P/T = 1. These large variations imply that Figure [1] cannot 
be used to estimate the errors in any particular case (except if the rms ratio is close to 
unity). Rather, it should be used as a general guide to the reliability of the experiment. Any 
individual planet detection must be analyzed based on the actual data and the parameters 
derived. 



5. Discussion 

The results derived here are primarily of interest in regard to future astrometric missions 
such as GAIA and SIM. The first point is that planets with periods that are even slightly 
longer than the mission cannot be reliably detected unless they are many times more massive 
than the nominal thresholds of detection. And second, if one is forced to allow for uniform 
acceleration, then the same statement applies to P/T > 0.7. Hence, to preserve the discovery 
space at the longest accessible periods, it is critically important to identify or rule out 
distant perturbers by supplementary data, such as longer-term RV observations to look for 
large, distant planetary or stellar companions, or perhaps AO observations to find stellar 
companions. This should be relatively straightforward for SIM, with its limited number of 
high-precision targets, but may be more difficult for GAIA, which is a survey instrument. 
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Fig. 1. — Panels (a) and (b): ratios of true amplitude and period error to naive eqs. (|Sj) and 
(jHJ) as a function of P/T and averaged over all phases <fi. Case of star with single planet 
is shown for RV {green curve) and astrometry [red curve). They differ because astrometry 
has one extra nuisance parameter (proper motion). If the star suffers an extra source of 
uniform acceleration from another companion, there is an additional nuisance parameter for 
both RV {red curve) and astrometry (blue curve). In all cases, RV follows the naive equation 
to higher P/T and degrades thereafter more gracefully. Panel (c): same ratio for the mass 
error relative to eq. flSJ). Both the isolated (dashed curves) and uniform-acceleration (solid 
curves) cases are shown for RV (blue curves) and astrometry (red curves). Note that the 
isolated-planet astrometry curve is no longer identical to the uniformly-accelerated RV curve, 
as was true for amplitude and period. This is because the mass is derived from different 
combinations of the amplitude and period for the astrometric and RV cases (eqs. [TDUTT] ). 
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Fig. 2. — Ratio of the true mass error to that given by the naive eq. (jHJ) as a function of 
phase for the special case P/T = 1. When the ordinate in Fig. [1] is large, the variations 
with phase are also large. Hence, Fig. Q] should be regarded as a guide to the statistical 
properties of the experiments rather than a precise estimate of the errors in any particular 
planet detection. 



